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The power spectrum of finite-temperature quantum electromagnetic fluctuations produced by 
elementary charge carriers under the influence of external electric field is investigated. It is found 
that under the combined action of the photon heat bath and the external field, the low-frequency 
asymptotic of the power spectrum is modified both qualitatively and quantitatively. The new term 
in the power spectrum is inversely proportional to frequency, but is odd with respect to it. It 
comes from the connected part of the correlation function, and is related to the temperature and 
external field corrections to the photon and charge carrier propagators. In application to the case 
of a biased conducting sample, this term gives rise to a contribution to the voltage power spectrum 
which is proportional to the absolute system temperature, the charge carrier mobility, the bias 
voltage squared, and a factor describing dependence of the noise intensity on the sample geometry. 
^"V It is verified that the derived expression is in agreement with the experimental data on 1/ /-noise 

measurements in metal films. It is shown also that the obtained result provides a natural resolution 
to the problem of divergence of the total noise power. 
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I. INTRODUCTION 

The origin of nicker noise [l| observed in virtually all conducting media remains an open issue in condensed matter 
\q , theory. Despite numerous models suggested since its discovery eighty years ago there is presently no consistent theory 
which would explain the main characteristic properties of this omnipresent noise. The power spectral density of 
""^5 \ flicker noise is proportional to J -7 , where / is frequency, and the exponent 7 is about unity (usually 7 takes values 
0.9 to 1.5). One of the essential difficulties for theoretical explanation is the fact that experiments show no limits for 
this dispersion law, neither lower or upper. Although flicker noise dominates only at sufficiently low frequencies, the 
l//-component is detected in the whole measured band up to 10 6 Hz. It is also experimentally established that the 
noise power spectrum of a sample is proportional to the applied bias squared, and roughly inversely proportional to 
Zj , its volume. 

Ubiquity of flicker noise and universality of its properties suggest existence of a simple reason for its occurrence. It 
\ is natural to expect this reason to have a quantum origin. Although some of the models suggested so far do consider 
various quantum effects as underlying mechanisms of flicker noise (such as, for instance, trapping of charge carriers), it 
may well be that its origin is to be sought at the most fundamental level. Namely, it is plausible that the phenomenon 
of flicker noise has its roots in the very quantum nature of interaction of elementary charges with electromagnetic 
field. From this point of view, the problem was attacked by Handel 0], who suggested that the observed flicker noise is 
related to the spectrum of low-energy photons emitted in any scattering process, which, according to Handel, has the 
1// profile and reflects the well-know property of bremsstrahlung, namely, the infrared divergence of the cross-section 
considered as a function of the energy loss. Later, the argument was modified and the so-called coherent quantum 1/ / 
effect described Q on the basis of quantum electrodynamic results of Kibble and Zwanziger 0,0 • Although Handel's 
theory was severely criticized in many respects [f| @], it has found support in independent investigations of Refs. 
Handel's approach is based on consideration of current fluctuations. An essentially different quantum approach to 
the problem was recently proposed [lol [r"0 | , in which flicker noise is treated as a voltage fluctuation originating from 
quantum fluctuations of individual electric fields of charge carriers. In this case, the noise power spectrum is found 
by evaluating the two-point correlation function of the Coulomb potential of an elementary particle, dispersion of 
this function being related to the quantum spreading of the particle wave packet. However, the magnitude of noise 
induced by an external electric field, given by this theory, turned out to be too small to explain the observed noise 
level. 
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The purpose of this paper is to show that the value of induced quantum noise is actually significantly higher than 
the previously calculated. It turns out that in evaluating the effect of external electric field it is essential to take into 
account statistical properties of electromagnetic field. The point is that simultaneous account of the effects of photon 
heat bath and external field leads to appearance of a contribution of new type. Namely, the new term in the noise 
power spectrum is odd with respect to frequency, and hence corresponds to the part of the correlation function, which 
is odd with respect to the difference of its time arguments. The underlying reason that makes the appearance of this 
term possible is the inhomogeneity in time of fluctuations produced by individual charge carriers. As a consequence of 
this inhomogeneity, the correlation function of each individual contribution to the Coulomb field fluctuation depends 
on both time arguments separately, rather than on their difference only, so that symmetry under interchanging of the 
arguments does not forbid appearance of the odd contribution. The time homogeneity is restored only after summing 
up all independent contributions. The new contribution stems from the connected part of the correlation function of 
particle's field, rather than from the disconnected one that was in focus of Ref. [ll|, and is related to the temperature 
and external field corrections to the photon and charge carrier propagators, in contrast to considerations of Ref. 
where thermal bath and external field affected only real particle states. 

The paper is organized as follows. In Sec. Ill Al we briefly discuss the role of the heat bath effects in evaluating the 
mean value and the correlation function of electromagnetic field produced by elementary particles, and identify the 
contributions relevant in the low-frequency regime. The power spectral functions of the Coulomb field and voltage 
fluctuations are defined in Sec. Ill Bl and written in the form convenient for explicit calculations which are carried out 
in Sec. IIIII The low-frequency asymptotic of the power spectrum of Coulomb potential fluctuations in the absence of 
external electric field is evaluated in Sec. IIII Al and is found to exhibit an inverse frequency dependence. However, 
the 1// part of the spectrum is cancelled in the expression for the voltage spectral density. The non- vanishing 
l//-contribution to the voltage power spectrum is obtained in Sec. IIII Bl upon account of the influence of external 
homogeneous electric field on the virtual charger particle propagation. Application of the obtained results to solids 
and comparison with experimental data is given in Sec. IIII Dl Gauge independence of our treatment of electromagnetic 
fluctuations is proved in the Appendix. 

II. PRELIMINARIES 
A. Heat bath contribution to particle propagators 

Consider a quantized system of charged particles interacting with electromagnetic field. Let T be the absolute 
temperature of the system. We are interested in the influence of finite temperature on quantum properties of the 
electromagnetic field produced by charged particles. Specifically, finite-temperature correlations in the values of the 
particle's Coulomb fields will be investigated. To the leading order in the electromagnetic coupling, these correlations 
are a single-particle effect, in the sense that in this case only fields produced by one and the same particle correlate. 
In what follows, we thus confine ourselves to systems which allow perturbative treatment of charge carrier collisions, 
either in terms of original particles, or in terms of quasi-particles (e.g., conduction electrons in metals). In the latter 
case, the particle's mass and energy-momentum relation should be replaced by the effective ones. 

In this section, we shall discuss some general features of the temperature effect, related to the heat bath influence 
on virtual propagation of the electromagnetic and charged field quanta. Evidently, the heat bath has no effect on 
the mean electromagnetic field of a charged particle. Indeed, under above assumptions, this field can be represented 
as the amplitude of one-photon emission in a transition between free charged particle states, contracted with the 
photon propagator. Since the 4-vector of momentum transfer to a free massive particle, p, is always spacelike, p 2 < 0, 
distribution of real photons appearing in the definition of the photon propagator is immaterial in the calculation of 
the field. As a consequence, quantities built from the mean field, such as the disconnected part of the correlation 
function, are not affected by the photon heat bath (details of evaluation of this part see in Refs. [l(| (HI)- Things 
change, however, when the connected part of the two-point correlation function of electric potential is considered. It 
is defined by the following symmetric expression 

C c Q ° n (x;x>) = ±(m\A (x)A Q (x') + A (x>)A (x)\in) , (1) 

where x and x' are the spacetime coordinates of two observation points, Ao is the scalar potential Heisenberg operator, 
and | in) denotes the given in state of the system "charged particle + electromagnetic field." In the two-photon 
processes, photon momenta are allowed to take on lightlikc directions, and hence the photon heat bath does contribute 
to the function C§g n (x; x'). 

As is well-known, the ordinary Feynman rules of the S- matrix theory are not generally applicable for the calculation 
of in-in expectation values, and must be modified, e.g., according to Schwinger and Keldysh [l3[. This complication 
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was overcome in by rewriting Eq. ([1} in the form 

C^(x; x 1 ) = Re(out\T{A (x)A (x')}\m) . ( 2 ) 

which allows the use of the S-matrix rules. This transformation uses equivalence of the onc-particlc in and out states 
(and Hermiticity of the electromagnetic field operator). It is applicable to the present case as well, despite the fact 
that now the in state is not one-particle, because we are not going to take into account scattering processes in the 
heat bath itself. This can also be shown directly at the diagrammatic level following the route of transformations 
taken in Ref. [14j . which is formally the same in zero- and nonzero-temperature cases. Thus, in order to calculate 
the in-in expectation value ([T]) taking into account the heat bath effect, the standard finite-temperature- field-theory 
techniques can be used [l5j . Below we employ a version of the real time formalism, developed in [16], which is 
especially convenient in actual calculations since the momentum space propagators in this formulation do not involve 
the step function. 

The real time formulation involves doubling of all fields, which will be specified by a two- valued lower index. 
According to the diagrammatic rules derived in [l6l ]. the photon propagator has the following matrix structure 1 



JVC*) = /(0i^(fe) e -^, ^(fe) = 4 7 r^(^;^ D ™W 



(3) 



where 



l(*0 D 2 2(k) 

D n (k) = -D* 2 {k) = - , D 12 (k) = D 21 (k) = j |fc J _ 1 , 

[3 = 1 /T being the inverse absolute temperature of the system. In applications to the problem of 1/ /-noise considered 
below, the value of the product /3|feo| turns out to be very small. For instance, even for frequencies as large as 
10 6 Hz, and temperatures as small as 1°K, it is less than 1CT 27 • 10 6 /(10 -16 ) = 1(T 5 (the factors 1(T 27 and lCT 16 are 
contributed by the Planck and Boltzmann constants, respectively), so the denominators in the above expressions can 
be replaced by /3|fco|, implying that the second term in Du dominates. On the other hand, the temperature effect 
on the propagation of massive particles is much less prominent. For instance, in the case of conduction electrons in a 
crystal (this case will be used throughout as a standard example), the particle energy e is of the order (h/d) 2 /m, where 
m is the particle mass, and d is the lattice spacing. Taking d ~ 10 _8 cm, m ~ 10 _27 g, we find that e/T ~ 10 +5 , so that 
the temperature contributions can be completely neglected (for fermions as well as for bosons), and the propagator 
taken in the simple diagonal form 

rf4fc v\<l>(u\„-ikx _ ( D fi(k) 



Df^k) = -Df;(k) = (m 2 -k 2 -iO)- 1 . 
This is for a scalar particle described by the action 



= d A x{{d^* +ieA^*){d^<P~ieA^(t))~m 2 ^^} , (5) 



where e is the particle charge. Account of particle spin, though adds some extra algebra, does not change the long- 
range properties of its field, so that following Ref. we work with the simplest case of zero-spin particles. The 
matrix propagators are multiplied in the interaction vertices, generated by the triple and higher order terms in the 
Lagrangian, with an additional minus sign for the product of 2-components, as in the Schwinger-Keldysh techniques. 
The arguments of the Green functions are treated as 1-component fields. Finally, external particle lines represent 
normalized particle amplitudes or their conjugates, according to whether the particle is incoming or outgoing, just 
like in the conventional techniques. 

As was mentioned above, the photon propagator is dominated by the temperature contribution (as long as the 
range of momentum integration contains lightlike directions, see discussion in Sec. MI A"| . while in the massive particle 
propagator this contribution is negligible. It is important, on the other hand, that the heat bath affects significantly the 
real particle propagation, i.e., external matter lines in the diagrams. Bilinears of the particle amplitudes representing 
these lines are expressed eventually via statistical distribution function (see Sec. IIII Bl for details). Thus, apart from 
explicit T-dependence coming from the photon propagator, the correlation function also depends on temperature 
implicitly through the particle statistical distribution. 



The photon propagator is written here in the Feynman gauge. The question of gauge independence of the correlation function is 
considered in the Appendix. 
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B. Power spectral densities of potential and voltage fluctuations 

Connected contribution to the power spectral density of electric potential fluctuations is obtained by Fourier trans- 
forming Eq. (TT]) with respect to the difference of the time instants t, t': 

+00 

C(x,x',t',uj) = J drC^ n (x,t' + T;x',t')e-^ T . (6) 

—00 

The upper and lower indices in the notation of the correlation function are suppressed in the left hand side, for brevity. 
We are interested ultimately in the power spectrum of voltage fluctuations, Cu, measured between two observation 
points x, x 1 (e.g., two leads attached to a conducting sample). The connected contribution to the voltage correlation 
function is given by 

Cu(x,x',t,t') = ~(m\U(t)U(t') + U(t')U(t)\m) , (7) 

where U(t) = Ao(x,t) — Ao(x',t) is the operator of voltage between the two points. This function is separately 
symmetric with respect to the interchanges x <-> x 1 , and t <-> t' , unlike the function CQQ n (x;x') which is only 
symmetric under x <-> x' . Substituting the definition of U(t) in Eq. ([7]), the former can be expressed via the latter 

Cu(x,x',t,t') = C$ n (x,t;x,t') + C$ n (x',t;x',t') 

- [C^ a (x, t; x', f ) + CggV, t; x, t')} . (8) 

Accordingly, the power spectral density of voltage fluctuations, defined by 

+00 

Cu(x,x',t',Lu)= J dTCu{x,t' + T,x\t')e-^ T , (9) 
—00 

is expressed through that of potential fluctuations as 

Cu(x, x', t', uj) = C{x, x, t', u) + C(x', x', t', uj) - [C(x, x', t', uj) + C(x\ x, t', uj)} . (10) 

Although Cjj(x,x' ,t,t') is symmetric with respect to the interchange t *-* t' , it depends on both time arguments 
separately, and therefore Cu{x 7 x' , t, uj) does not have to be an even function of uj. 

When calculating the power spectrum of potential fluctuations according to Eqs. @, ©, it is convenient to perform 
the Fourier transformation under the sign "Re" in Eq. For this purpose we introduce the Fourier transform of 
the two-point Green function: 

+00 

G(x,x',t',oj)= J drG{x,t' + T;x',t')e- lUT , G(x; x') = {out\T{A (x)A (x')}\m) , 

— OO 

(11) 

with the help of which the power spectral density of potential fluctuations can be written as 

C{x,x',t',uj) = ^Re{G{x,x',t' 1 uj) + G(x,x',t',-uj)} 

+ ^Jm{G(x,x',t',uj)-G(x,x',t',-uj)} . (12) 

We see that contributions to the function C(x,x' ,t' ,uj), and hence to the voltage power spectrum, are either real 
even, or imaginary odd functions of frequency. 

III. EVALUATION OF LOW-FREQUENCY ASYMPTOTIC OF SPECTRAL DENSITY 
A. Power spectrum in the absence of external electric field 



The tree contribution to the right hand side of E q. Q when the influence of external electric field is neglected 
is shown in Fig. [TJ Repeating the argument of Ref. [llj , it is not difficult to show that the leading contribution is 
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contained in the diagrams QJa) , [T^b) ■ 
have 



Distinguishing their contributions by the corresponding Latin subscript, we 



G a (x,x') = ie 2 J J d 4 zd 4 z' {D(x lZ ) 
G b (x,x') = G a (x',x), 



0o(z)5 o ^(z,z')9 o '05(z') 



D(z',x' 



(13) 
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where 



and 4>o is the given particle state. Going over to momentum space with the help of Eqs. ((3]), ((4]), introducing the 
spectral function for G a (x,x') according to Eq. (|11|) . and writing the matrix product longhand yields 



P° = £q+ P - e q , e q = +\Zq 2 + m 2 , 



(14) 



where 



J a (p,q,x - x',oj) = -i 



d 3 k 
(2tt) 3 



■ik(x' — x) 



(2q + k )(2q + k + p ) 



D n {k)Dl x (q + k)D u {k-p) + D 12 {k)D$ 2 (q + k)D 21 (k - p) 



kn—UJ 



Here is the particle 4-momentum, a(q) its momentum wave function at some time instant to •> normalized by 

d 3 q 



1 



K«)r = i 



(15) 



(16) 



and it is taken into account that Df 2 = . 

The second term in the square brackets in Eq. (JT5J) can be neglected. Indeed, in view of the factor -Di2(fc) which is 
proportional to S(k 2 ), and the condition fco = u>, the momentum k contributes only a tiny value to the argument of the 
factor D21; for electrons in a crystal, for instance, the ratio |fe|/|p| is of the order (fko / c) / (H / d) = lodjc ~ ui ■ 10~ 18 s. 
Therefore, this factor can be written simply as D-2\{p) ~ S(p 2 ). On the other hand, momentum transfer p to the 
massive particle is spacelike, and hence the argument of the delta-function is always nonzero. Furthermore, using 
explicit expression for the photon propagator, the first term in the square brackets reads 



1 2iriS(k 2 
k 2 + i0 ~ e^lM - 
27nS(k 2 ) n0 



^11 
1 



1 



{k-pf 



2niS((k-p) 2 

iO e P\k -po\ _ 1 

^2 



1 



1 ^11 



1 



[k 2 + iO] 11 [{k - p) 2 + tO] 
4 2mS((k - p) 2 ) 2mS(k 2 ) n0 2md((k - p) 2 ) 



+ 



0|fcol - 1 ^{k-p) 2 +iO k 2 +i0 11 e /3|fco-po| _ 1 e P\k \ _ I 



o/3|fco-pol — 1 



As before, the last term on the right hand side can be neglected, while the first term, describing the zero-temperature 
contribution, has already been considered in Ref. [111 ]. Furthermore, fco enters the temperature exponent in the third 
term in the combination (fc — Po), and, therefore, this term does not contribute to the leading order in l/u), because 
in practice \oj\ <§C po- Indeed, estimating the energy transfer as po ~ (pq)/m, and taking \p\ ~ l/l, where I is the 
characteristic sample length, one finds for our standard example |w/po| ~ 10 _8 |w|/, where l,u> are supposed to be 
expressed in the CGS system of units. Even for I as large as 1 cm, this ratio is very small for all practically relevant 
frequencies. Thus, the contribution of the second term only remains to be considered. In view of the factor <5(fc 2 ), the 
pole of the function Dn(q + k) in this term does not contribute. This is again a consequence of the requirement that 



the momentum transfer to a massive particle on-shell be spacelike: conditions fc 2 = 0, q 2 = m 2 , and (q 
cannot be satisfied altogether. Hence, the scalar particle propagator can be written simply as Dn(q + fc) = 
or, assuming that the particle is nonrelativistic, \q\ <C m, as 



fc) 2 



DUq + k) 



1 

2muj 



-l/(2qk), 



(17) 
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Finally, neglecting k in comparison with p in the factor l/[(k — p) 2 + iO], using p 2 < to omit iO, and retaining only 
terms singular in u>, we find 

uj(3\uj\p z J {2iry 
(2m + p ) 



2ttujPp 2 



0(u) . (18) 



The reason why po has been kept along with m in the numerator will become clear soon. As to the contribution of 
the diagram [IJb) , changing k — » fc + p in Eq. (fT5|) , and then q — > g p, p — > p in Eq. (fT?)) shows that 

G b (x,x',lf,u) = G* a (x,x',t',-cj). 

Thus, the total contribution to the function G(x, x' , t', u>) is 

G(x,x',t',Lo) = G a {x,x',if,u)) + G* a { X , X ',if,-w). (19) 

It is seen from this relation and Eq. lfl8|) that only imaginary part of G a gives rise to a nonzero contribution to the 
total Green function G, and this part corresponds to the term proportional to po on the right hand side of Eq. (|18[) . 
Thus, 

G{x,x>,t>,») = [[ ^^ fl M«*talp)^e«»(-»)-', (20) 

m^pui J J (27r) J (2tt) 6 p z 

where, in the denominator, the energies e q ,e q + p have been replaced by m, and p° neglected in comparison with |p| 
on account of the condition \q\ <C m, while in the numerator, p° has been replaced by its leading long-range term, 
(pq)/m. It is instructive to see what the right hand side of Eq. (|20p becomes in a particulary simple model case when 
the amplitude a(q) can be written as 

a(q) = b(q)e-^> , (21) 

where b{q) is a real function of the particle momentum. In this case, x$ is easily identified as the mean particle 
position, and hence, the function 6(g) describes the momentum space profile of the particle wave packet. After 
extraction of the position-dependent phase factor, the amplitude becomes a relatively slowly varying function of the 
particle momentum, therefore, to the leading order of the long-range expansion, 6(g +p) in the integrand of Eq. (j2"D|) 
can be replaced by 6(g) 

G(x,x',t',u) = -i^- // f« f P ^(g^Me^'-toH^-a,) , (22) 
m z puj J J {Ziry (27r) J p* 

In the exponent, po generally cannot be replaced by (pq)/m, because the subleading term of its long-range expansion, 
p 2 1 (2m), although small compared to (pq)/m, can change the phase significantly, provided that the difference (t' —to) 
is sufficiently large, as is the case if the particle collisions are neglected. However, taking into account the latter reduces 
this difference to the particle mean free time, r/, so that the product po(t' — to) can be completely neglected. Indeed, 
Tj can be estimated roughly as d/(\q\/m), and hence, po(t' — to) ~ |p||g|r//m ~ d\p\ ~ d/\x' — x \ <C 1. Then the p 
integral can be evaluated with the help of the formula 



d3 P _i<j*)*"P 



IX 



3 «2 |™|3 



yielding 



(2tt) 3 p 



G(x,x',t',uj) = -2t^- { -^, r' = x'-x , (23) 



where the overline denotes g-averaging over the given particle state. In the absence of external electric field, q is zero, 
and therefore so is the right hand side of Eq. (f2"3")l . As we will see in the next section, the same result is obtained in the 
general case without the use of the model decomposition (f2"Tj) . But even for nonzero g, the voltage power spectrum 
calculated from G given by Eq. ([23)1 turns out to be zero. This is verified directly by substituting expression (|23|) into 
Eqs. (fT2"|) , (|10p . The reason for nullification of the voltage power spectrum is easily identified - it is the consequence 
of the fact that the function G a (x, x', t', u>) is independent of the x coordinate. The ^-dependence has been lost upon 
extracting the low-frequency asymptotic of J a in Eq. JT| 
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B. The influence of external electric field and particle collisions 

Let us now consider corrections to the power spectrum due to constant homogeneous external electric field, taking 
into account also the influence of particle collisions. These corrections are twofold. First of all, the field affects the 
particle wave function a(q), which is symbolized in Fig.[^]by inserting the vertices of particle-field interaction into the 
two external solid lines. This is only a schematic picture, because the effect of a constant homogeneous field on the 
free particle states cannot be treated perturbatively. The latter circumstance, however, is not important in view of the 
particle collisions which prevent the particle from gaining too much momentum from the field, thus cutting down its 
effect (particle collisions are symbolized in Fig.[5]by a virtual photon interchange between particles). Account of these 
two factors is accomplished by replacing the particle momentum probability distribution, a(q>)| 2 , by the statistical 
distribution function, obtained as a solution of the kinetic equation in the presence of external electric field. This 
point will be discussed in more detail later in this section. 

Next, the particle propagator is also modified by the external field. It is not difficult to see that for a sufficiently 
small field strength, E, this modification can be treated perturbatively. It can be recalled that in coordinate space, 
it amounts to multiplying the zero-field propagator D^{z, z') by 2 exp {ie(E 7 z + z')(z' — z )/2} . Although the low- 
frequency limit is determined by the large-time behavior of the quantities involved, implying that (z' — z Q ) ~ 1/lu is 
large, the exponent can be made as small as desired for any given u> by taking \E\ sufficiently small. This also will be 
clear from the explicit calculations to follow. The lowest order correction to the correlation function is represented by 
diagrams with a single insertion of the particle-field interaction vertex into the internal solid line, as shown in Fig. [3] 
As in the preceding section, these diagrams vanish unless all interaction vertices are 1-type, on account of momentum 
conservation in the vertices together with the mass shell conditions for the massive particle. The contribution of the 
diagram [3ja) to the spectral density of the two-point Green function 



G*(x,x',t',Lo)= (47re) 



d 3 q d 3 p d 3 ki a(q)a*(q + p) 



(2tt)3 (2tt)3 (2n) 3 ^/2s q 2e q+p 
xe v ° {t '- to] - l(p - kl)x ' ipik^J^ {p, q,x- x', u) , (24) 



where 



/d 3 k i 
__ eJ fc(* - x){2qo + ko){2qo + 2ko )(2q + k Q +po) 

D ll {k)Di l {q + k)Di 1 {q + k + k l )D ll {k + k 1 -p) , fc 1 = (0,fci), (25) 

ko—uj 

and <p(k) is the Fourier transform of the external field potential, <p(x) = —(E,x). We do not include an arbitrary 
constant in this expression because it is clear in advance that it cannot affect the final result. It is proved in 
the Appendix that the correlation function is actually invariant under the most general gauge variations of the 
electromagnetic potential. Substituting 

^(fc) = -*(2tt) 3 ( E -^S^(k) 

into Eq. I|24[) and integrating by parts, the k\ integral in Eq. is brought to the form 

i J ^k^Qc) [Dt 1 (q + k + k 1 )D 11 (k + k 1 ~p)e- l ^- k ^ x '] , (26) 

where only terms involving k± are retained. As we have seen, the singularity at u> = in the function G a (x, a;', t', uj) 
comes from integration over small k, and this singularity is now strengthened by the extra particle propagator and 
the differentiation with respect to fci. 3 Consider first the case when the fci-derivative acts on the last two factors in 
the square brackets. Since these depend on the difference (p — k±), changing d/dki — > —d/dp, and then integrating 



2 This phase factor incorporates contributions involving only vertices linear in the photon field. This is sufficient for the subsequent 
discussion concerned with the leading correction to the propagator. Inclusion of the other interaction vertex adds terms of higher orders 
in E to the exponent. 

3 Independently of the proof given in the Appendix, it can be noticed that the constant term in the potential does not involve the 
fci-differentiation, and hence does not contribute to the leading singularity for ui — > anyway. 
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by parts with respect to p in Eq. ([24]) . this derivative is rendered to act on terms independent of fc 4 . Thus, of the 
three factors in the square brackets, only the first is to be differentiated in effect, thus reducing the expression to 



2i(E,q + k) 

[m 2 — (q + fc)]' 



D n {k - p)e-^ px "> 



Taking into account also the vertex factor e(2q + 2fc ), we see that the first order correction to the charged particle 
propagator due to constant homogeneous external electric field is obtained by inserting the factor 

iie(q + k )(E,q + k) _ 
[m 2 — (q + fc)] 2 

into the integrand in Eq. (p~5|) . The rest of the calculation repeats the steps of Sec. MI Al Substituting explicit 
expressions for the photon propagators, one sees that only the part proportional to 5(k 2 ) is to be retained in the 
expression for Du(k), while the corresponding part in Dn(k-p) is to be omitted. One consequence of this observation 
is that the factor k simplifies to 

ie(Eq) 

mto 2 

where we have taken into account that fco = u> <C qo, and |fe| = \ui\ <C \q\ [indeed, taking our standard example of 
electron in a crystal, the ratio |w|/|q| is, in the ordinary units, (7l|cj|/c)/|<j| ~ \u)\d/c ~ |cj|10 -18 , with u> expressed in 
Hz]. Another consequence is that the contributions of diagrams [3ta), GUb) are related, as before, by 5 

G?(x,x',1?,w) = G**(x,x',1/,-l>). 



Furthermore, Eq. (|18p is now replaced by 



2ir(2m+p Q ) f _fk_ ik^'-x] 



uj(3\lo\p 2 
ie(Eq) (2m + p ) 
mto 2 2nuj(3p 2 



(2tt) 3 



5{lo 2 



1 



6 



{x' 



O{to) 



(27) 



An important difference in comparison with the result of the preceding section is that because of the extra imaginary 
unit brought in by the factor tt, the term proportional to po in the last formula gives rise now to a purely real 
contribution upon substituting in Eq. (|2^|) . and hence is cancelled 6 in the sum of diagrams Oa), [3Jb) . On the 
other hand, the terms independent of po survive, and lead to the following expression for the first order correction to 
the spectral density of the two-point Green function 



G E (x,x',t',uj) 



1 



xRe 



d3q ^(g+pjW^w 



(2tt) 3 (2tt) 3 



P 



(28) 



In a many-particle system, this result is to be expressed through the one-particle density matrix, which is accomplished 
by replacing a*(q')a(q) — > Qo(q', q), where go is the momentum space density matrix at the time instant in- We recall 



4 Moreover, if the function a(q) is decomposed as in Eq. i'2\\) . and the potential is chosen respectively to vanish at the point tco, i.e., 
ip(x) = —(E, x — xq), then the p-derivative acts on the product of e IpK ci with the slowly varying factors b(q+p), £ q + p , po- Differentiation 
of the exponent gives xo which just cancels the purposely chosen constant term in the potential, while the result of differentiation of 
the remaining factors can be neglected to the leading order of the long-range expansion. This observation can be useful in assessing the 
higher order corrections to the correlation function. 

5 It is convenient to prove this relation before the integration over fci using the following sequence of substitutions: fc — > fc + p — fci , 
q —* q — p, and then p — * —p. The extra factor (—1) coming from complex conjugation of the imaginary unit in the factor >c 
is compensated by that from the integration by parts with respect to fci. Alternatively, this can be proved directly in coordinate 
space (after extracting the relevant contribution) by substituting exp \ie(E, z + z')(z' — zo)/2} D^{z, z') for the particle propagator, 
interchanging the integration variables z, z' , and taking into account symmetry of the functions D(x, z), D'f'(z, z') . 

6 Here p is neglected in comparison with q. Otherwise, there is a residual term proportional to po in the total of the two diagrams. The 
usual manipulation with the integration variables q q p, p — » — p shows that this term can be obtained by substituting (Eq) — > (Ep). 
Its ratio to the main contribution considered in the text is |p|po/™|<z| ~ p 2 /m 2 . For electrons in a sample of characteristic size icm, 
this is, in the ordinary units, (h/lmc) 2 ~ 10 -20 /I 2 . 
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that to denotes the instant at which the particle state a(q) is prepared. It can be identified, for instance, as the moment 
the charge carrier enters the sample, or escapes from a surface trap, etc. The factor e lpo ^ ~ to ' in the integrand then 
realizes time evolution of the density matrix from to to t'. Since po = (q + p) 2 /2m — q 2 /2m, this is a free evolution. 
The fact that the evolution of charge carriers in solids is not actually free on macroscopic scales is not important for 
the present consideration in which appearance of the 1 /w-singularity is related to the effects of the medium on the 
field propagators. In this respect, it is essentially different from considerations of Ref. [ll[, in which time evolution 
of the particle wave packet was the central issue, and, in particular, the requirement that the particle collisions be 
elastic was important. To take into account particle collisions in the present case, it is sufficient to consider them as 
instantaneous. Then the interval (to, t 1 ) is divided into a sequence of short time intervals of order ry (the particle mean 
free time), on each of which the density matrix evolves freely as in Eq. |28|) . and changes abruptly at the collision 
instants. Going through this sequence the density matrix tends to the stationary statistical distribution function, 
q(q' ' i q)> which is independent of the initial particle state. What is important here is the sign of the difference (t' — to). 
Recall that t' is a fixed time instant to count off the time interval r with respect to which the correlation function 
is Fourier-transformed, and that each particle has its own to. This means that for a given u>, the system is observed 
during the time interval (t' — At,t' + At), where At ~ 1/u, and tos are distributed uniformly over this interval. The 
density matrix evolves forward (backward) in time, if t' > to (t' < to). But time reversal involves inversion of particle 
momentum, and therefore, the reciprocal contributions to the function G E (x, x', t 1 , uj) have opposite signs. To be 
more specific, let {t' — to) > 0. Then the exponent e lpo ^ ~ to > realizes a forward evolution of the density matrix, so 
that the integral in Eq. ([2"5)l takes eventually the form 

^ d3P : g (q + P ,q) { -^e^'K (29) 



(2tt) 3 (2tt) 3 ^ p 

On the other hand, if (t' — to) < 0, then the density matrix evolves backward. In momentum space, the initial state 
of the reversed motion is represented by the amplitude a(q) = a*(—q). Taking complex conjugate of the integral in 
Eq. ([28]) (which does not change the value of G E in view of the sign "Re"), and changing the integration variables 
q — > q, p — > p gives in this case 

d3q d3p r a( g )a*(g+p)^e^(«o-o-i(i«') . 



(2tt) 3 (2tt) 3 v ^ " p 

After replacement a*(q+p)a(q) — > g(q+p, q), where g plays the role of the momentum density matrix at the moment 
t' , the exponent e m ' t - t °~ t ' governs forward evolution of this state on the interval (t' , to), so that the above expression 
takes the form 



(2tt) 3 (2tt) 3 ^ p 

The density matrix here is the same as in ([2"9")l . because the statistical distribution is independent of the initial state. 
We see that reciprocal contributions to the function G E (x, x 1 , t', uj) cancel each other when summed over all particles 
in the system. Thus, we arrive at the important conclusion that the total noise intensity is independent of the number 
of particles, and remains at the level of individual contribution. As was shown in Ref. [ll| , this conclusion is also true 
of the disconnected part of the correlation function, though by virtue of quite different reasons. 

It is customary to further express the function g(q',q) via the real mixed distribution function, n(r,q), according 

to 



g{q +P,q) = J d 3 re l 



(pr) n (r,q+^ 



Probability distributions for the particle position in a sample and for its momentum can be obtained by integrating 
n(r,p) over all p and the sample volume, respectively. Using this in the expression ([29]) , and substituting the latter 
into Eq. (QSJ) yields 

G E {x,x',t',u J = — 1-— (x'-x) 2 

-,///^^<^„(„ + f),— . (30, 

As we know, the first term in the square brackets in this formula doest not contribute to the voltage power spectrum, 
and can be omitted. Furthermore, after shifting q — > q — p/2, and omitting the imaginary term proportional to (E,p), 
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the triple integral becomes purely real, so the symbol "Re" can be omitted. Integrating then over p with the help of 
the formula 

d3 P _ e i(P"0±[ = _L 



(2tt) 3 p 2 \x\ 



G E (x,x\t',u,)= 2ie y- x)2 ff 4^dMEq)^4j. (31) 



we thus obtain 

H E (cr. r> t' = — 

3m/3uj J J (27r) 3 \r — x' 

Equation (Ti"2"]) shows that the spectral density of the correlation function is given by the same expression (|3"Tj) , so that 
the power spectrum of voltage fluctuations is, according to Eq. (TIT)]) . 

- ~ 2J %^ I ^mr)) ( ' + , (32) 



where 



is the local drift velocity of the charge carriers, f2 denoting the sample volume. For a crystal in the homogeneous 
external field, v is a function of the crystalline direction, 

Ui = [iikE k , i,k = 1,2, 3, 

where /i^. is the charge carrier mobility tensor. With the help of this formula Eq. (f32|) can be rewritten as 

C u (x,x>,t',u;) = -ir l l £, r, = ^-, (33) 
uj 3p 

where 

E 

Uq = \E\\x — x'\ is the bias applied to the sample (it is assumed that E || (x — x'), as is usually the case), and g is a 
geometrical factor 

^U d3r (\^\ + vh7])- (34) 

If Fourier transformation is defined in a purely real form, i.e., as a decomposition in cos(ojt), sin(ajT), rather than in 
e luJT , then the spectral density is also real: 

C v {x,x',t,u,)=r 1 Vl, V = ^-, (35) 
uj 3p 

We mention for future reference that if the sample is an elongated parallelepiped with the leads attached to its ends 
(as is usually the case in practice), then the g-factor can be evaluated as 

2 f whdx 2 I 
g « — / = - In - , 36) 

Iwn J X i w 

w 

where I, w, h denote the sample length, width and thickness, respectively, and it is assumed that h < w <C /■ We note 
also that in the ordinary units, the dimensionless factor ry reads 

2e 3 fig 2a 2 

a being the fine structure constant. In particular, in the CGS system of units, 

r, w 3.4- 10- 22 c^T, (37) 
where the absolute temperature T is to be expressed in °K. 
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C. On the unboundedness of l//-spectrum 

In this section, a special feature of the derived expression for the power spectrum, namely, its oddness in frequency, 
will be discussed in connection with the problem of observed absence of frequency limits of the l//-law. As was 
mentioned in Introduction, 1// noise has been detected in a very wide frequency band ~ 10 _6 Hz to 10 6 Hz. This 
fact represents one of the essential difficulties for theoretical explanation, because all physical mechanisms underlying 
existing models of flicker noise work in much narrower subbands, and none of the models suggested so far has been 
able to explain the observed plenum of the l//-spectrum. 

On the other hand, existence of bounds on this spectrum is generally believed to be necessary in order to guarantee 
finiteness of the total noise power. There is a well-known argument [171 ] according to which these limits are actually 
unnecessary when the flicker noise exponent 7 is strictly equal to unity, because the logarithmic divergence of the 
total power is not a problem in this case in view of the existence of natural frequency cutoffs such as the inverse 
Planck time and lifetime of Universe. However, this reasoning does not work for 7 ^ 1, in which case divergence is a 
power of the cutoff. At the same time, the results obtained above reconcile unboundedness of 1/ /-spectrum with the 
requirements of stationarity and finiteness of the total noise power in a quite natural way. Indeed, using Eq. (|33p we 
find 

— OO 

More generally, if the spectrum C\j ~ l// 7 is continued to negative /'s as an odd function, then for any < 7 < 2 
the integral 

+00 +00 

-00 

is convergent in both limits / — > and / — > oo. In particular, the singular contribution to the voltage variance (i.e. 
to the quantity Cu(t, i')| T=0 ) vanishes. 

Since appearance of odd contributions to the power spectrum is somewhat unusual in macroscopic fluctuation theory, 
let us discuss it in more detail. Under stationary external conditions, the voltage noise power spectrum (to be denoted 
below simply as Cu(t,t'), with the spatial arguments suppressed, for brevity) must be independent of t' . This is an 
expression of the noise stationarity, or, using a term more suitable for the subsequent discussion, time homogeneity 
with respect to the macroscopic system. It is usually realized as the requirement that Cjj(t,t') be a function of 
the difference t — t' = r. Since Cu(t,t') is also symmetric with respect to the interchange t <-> t', an immediate 
consequence of this is that it is actually a function of |t|, and hence the spectral density is a real even function 
of frequency. It is important, on the other hand, that time homogeneity is not necessarily exhibited by individual 
contributions to the total voltage fluctuation, whatever mechanism of flicker noise generation be. In particular, this 
property evidently does not take place at the microscopic level, i.e., with respect to elementary processes such as 
charge carrier trapping, surface or grain boundary scattering, etc. Stationarity of the macroscopic process emerges 
usually upon summation over a large number of individual contributions, so that this microscopic inhomogeneity turns 
out to be inconsequential. However, this summation is not the only way to obtain a stationary correlation function 
symmetric in t, t 1 . Another possibility, which is realized in the present paper, is that flicker noise may be a one-particle 
phenomenon, in the sense that the entire effect can be ascribed to elementary fluctuations produced by single charge 
carriers. In this case the function Cu(t,t') does not have to depend solely on |r|, and as the explicit calculations of 
Sec. IIIII show, it actually does not. As was mentioned above, elementary processes are inhomogeneous in time, and 
hence the symmetry with respect to t *-+ t' imposes no restriction on the r-dependence of the correlation function. 
The only remaining requirement, namely reality of the correlation function, implies that contributions to the spectral 
density must be real even, or imaginary odd functions of frequency [Cf. Eq. (|12| ], These two cases correspond to 
the Fourier decomposition of the function Cu(t' + r,t') in cos(wr) and sin(wr), respectively, and describe the parts 
symmetric and antisymmetric with respect to the difference of its time arguments. Finally, transition to the statistical 
distribution removes the ^-dependence of the power spectrum [Cf. transition from Eq. I|28p to Eq. ([3H)) ] . This restores 
macroscopic time homogeneity of the correlation function, but leaves the possibility of being odd with respect to 
the difference of its time arguments. In other words, dependence of the power spectrum on t' shows itself only at 
microscopic scales, while macroscopically fluctuations look as if they were homogeneous in time. 

The l//-spectrum derived in the previous section has no lower frequency cutoff. As to the upper bound, it is given 
by the condition / CT [see Sec. Ill Aj . or in the ordinary units, / <C kT/h » 10 11 THz, with T expressed in °K. We 
see that from the practical point of view, the obtained spectrum has no upper cutoff either. 



dui sin(o;T) r r]U( 
2^ uj ~ H ~T 



12 



D. Validity of Eq. (135[) and comparison with experimental data 

Let us next discuss the range of applicability of the obtained results. The validity of the perturbative treatment of 
the external field imposes a very strong bound on the field strength. For this purpose we first collect all characteristic 
factors that have appeared in the course of extracting the 1/w-asymptotic of the voltage power spectrum. As we 
have seen, insertion of the vertex describing interaction of the virtual charged particle with external field amounts to 
multiplying the zero-field diagram by the factor, in the ordinary units, e(Eq) / '(muj 2 h) , which is eventually promoted 
by the q-integration into e^,E 2 / (u> 2 K). Furthermore, the leading non- vanishing contribution to the voltage correlation 
function has been obtained after expanding e lk ( x - x ) i n the integrand of the fc-integral, which brought in a factor 
u> 2 (x — x') 2 /c 2 . Thus, the overall factor is efiE 2 (x — x') 2 /(c 2 h) — efiUQ /(c 2 h). This is small provided /iUq <C c 2 h/e w 
10 3 units CGS, which is a quite soft requirement met in virtually all flicker noise measurements (Cf. examples below). 
The problem with this estimation, however, is that at higher orders, the scalar product (Eq) is to be estimated 
as |-E||<z|, because (Eq) 2 is of the order (|i£||<7|) 2 , rather than (Eq) . As a result, the requirement that the factor 
e(Eq)/(muj 2 h) be small leads to the following upper bound on the electric field strength for a given frequency w, in 
the CGS system of units, \E\ <C draio 2 /e w 10~ 25 cl> 2 . At the same time, the values \E\ ~ 1 are quite normal in flicker 
noise measurements. In other words, from the point of view of the developed theory, the experimentally relevant 
regime is identified as the strong field limit. Yet the use of Eq. ([55]) in this limit can be justified to a certain extent 
by recalling that the perturbative expansion is in reality an asymptotic expansion, and hence the fact that Eq. (|35[) 
gives the first non-vanishing term of the voltage noise power spectrum implies that the question of validity of the 
perturbative expansion is actually a question of whether or not it is legitimate to use this expansion to obtain higher 
order corrections to Eq. (|35[) . A rigorous justification is a difficult task because it requires the use of non-perturbative 
methods. Thus, this issue is left open until careful investigation of the strong field limit. One of the possible ways 
this problem can hopefully be resolved is a partial summation of the perturbation series, followed by an analytical 
continuation with respect to k. 

After this discouraging observation of strong divergence of the asymptotic series, the more striking turns out to be 
the fact that Eq. (|35p is in a general agreement, qualitative and even quantitative, with the existing experimental data. 
First of all, the spectral density is quadratic in the applied bias. This is perhaps the most solidly established property 
of flicker noise. Second, the noise level is inversely proportional to the sample size. As to the dependence of flicker noise 
amplitude on sample dimensions, agreement in the literature is not that good. Experiments are usually arranged so as 
to prove one of the two main competing points of view on the flicker noise origin, namely wether it is a bulk or surface 
effect. Although this issue is far from being resolved, there is no doubt that the noise level increases with decreasing 
sample size. Third, it is generally agreed that, with other things being equal, the flicker noise is more intensive 
in semiconductors than in metals, and this is again in conformity with Eq. (|35[) . because charge carrier mobility 
is higher in semiconductors than in metals, usually by several orders. Unfortunately, determination of mobility in 
semiconductors (or semimetals) is a difficult problem, both theoretically and experimentally, and different experiments 
often give significantly different results. By this reason, the subsequent quantitative consideration will be carried out 
for metals only. Even in this case careful estimation of the noise level takes some effort. This is because electron 
mobilities in thin metal films commonly used in flicker noise measurements differ essentially from the corresponding 
bulk values, varying non-monotonically with the film thickness, and exhibiting complicated temperature dependence. 
Thus, the thicker the film, the more reliable comparison of theoretical and experimental results. Fortunately, the 
modern instrumentation allows measurements in sufficiently thick samples, electrical transport in which has bulk 
properties (usually, effects related to film thickness become important for h less than a few hundred nanometers) . As 
is well known, temperature dependence of the electron mobility in this case is well approximated by the 1/T law. 
Theoretically, this approximation is valid for T higher than the Debay characteristic temperature, but in most cases 
it is practically applicable already for T > 50°K. Thus, it follows from Eq. (|35|) that the flicker noise level in thick 
samples is temperature independent. This conclusion is confirmed, e.g., by the results of Ref. [l8T ] where 1// noise 
was measured in 2.44 /im thick metal films, which is quite sufficient for bulk treatment of the sample conduction. 
According to Fig. 5 of Ref. |l8|, the flicker noise level is constant for T > 50°K indeed. Unfortunately, the authors of 
[HI did not specify the metals used in their experiments, which makes further comparison with Eq. (|33|) impossible. 

In order to compare the absolute value of the noise spectral density given by Eq. (|35p with experimental data, 
we use the results of the classic paper [l!| where flicker noise in thin metal films was investigated. The information 
provided in this paper is sufficient for estimation of the noise intensity in the gold film shown in Fig. 2 of fl9j | . This 
was an elongated sample with h — 25 nm, w — 8 fim, I — 625 /im, biased at Uq = 0.81V, and operated at about 
40°K above room temperature. Substituting the sample dimensions in Eq. (|36[) gives g — 140 cm -1 . Estimation of the 
electron mobility is more subtle. As was mentioned above, charge carrier mobility in thin films strongly deviates from 
its bulk value, and this deviation is the main source of uncertainty in evaluating the noise level. In the case under 
consideration, fi is isotropic and can be found using the relation /i = a /en, where a is the electrical conductivity of 
gold, and n = 5.9 • 10 22 cm -3 is the free electron concentration. The bulk conductivity of gold at T = 330°K is equal 
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to 4.0 • 10 7 f2 _1 m _1 , but in thin films the value of a is strongly affected by the grain boundary and surface scattering, 
surface roughness and other factors. The relevant value of conductivity can be calculated indirectly using the I — V 
characteristic of the given gold sample, shown in Fig. 3 of 19]. According to this figure, the sample resistance was 
about 100 Q. Taking into account the sample dimensions given above, this implies that a = 1.2 • 10 6 r2 _1 m _1 . It should 
be mentioned that this value is approximately six times lower than that obtained in more recent studies of electrical 
transport in thin films. For instance, according to Ref. [22j conductivity of a 25 nm thick, 15 /j,m wide gold film obtained 
by a laser-improved deposition of nanoparticle suspension, is 7.1 • 10 6 f2 _1 m _1 . The same value can be obtained also 
indirectly using the data given in Refs. [13, Hl|. According to [2(|, the conductivity of gold is 75% to 85% of its bulk 
value for h = 100 nm, depending on the choice of the substrate, and decreases below that value approximately linearly 
with decreasing thickness. On the other hand, according to Ref. [2l[ conductivity drops to about 3 • 10 5 _1 m _1 for 
h = 5 nm. One readily finds from this that for h = 25 nm, a = (6.5 -j- 7.5) • 10 6 f2 _1 m~ 1 . Presumably, this difference in 
the values of conductivity is to be attributed to the quality of film deposition. In the case of a = 1.2 • 10 6 f2 _1 m _1 , 
the electron mobility equals to \i = 1.3cm 2 /Vs, and then Eq. (j37|) gives T] = 6.0 • 10~ 15 . Substituting this together 
with the bias value given above in Eq. (f35|) . and setting ui = 2ttJ yields Cjj = 6.3 • 10~ 16 V 2 /Hz for the frequency 
/ = 1 Hz, which is to be compared with the experimental value Cjj « 10~ 15 V 2 /Hz . 

IV. DISCUSSION AND CONCLUSIONS 

We have shown that the combined action of the temperature and external field effects results in appearance of a 
principally new contribution to the power spectral density of quantum electromagnetic fluctuations, given by Eqs. (I32p . 
(1331) . The power spectrum is thus modified both qualitatively and quantitatively. Being odd with respect to frequency, 
the new term in the power spectrum describes correlations in the values of voltage measured at two time instants, 
which are finite for all times. The underlying reason that makes the appearance of the new term possible (apart from 
the two factors mentioned in the beginning of this paragraph) is the inhomogeneity in time of fluctuations produced 
by individual charge carriers. As discussed in Sec. IIII CI oddness of the found l//-contribution gives a natural 
explanation to the observed unboundedness of flicker noise spectrum. Although the obtained result is valid, strictly 
speaking, only for very weak fields, we have seen in Sec. IIII 01 that it is in qualitative and quantitative agreement with 
experimental data even beyond its formal range of applicability. 

Next, an important qualitative difference of the present considerations from those of Ref. [Io| is to be emphasized. 
As we have seen in Sec. frequency dependence of the power spectral function is determined completely by internal 
structure of the Feynman diagrams representing the connected part of the correlation function. In other words, 
dispersion of the correlation function, considered in the present paper, is related to the properties of virtual quanta 
propagation, and not to the time evolution of the charge carrier wave function. This is in contrast to considerations 
of Ref. [l(| where the 1/ / asymptotic of the power spectrum was related to the spreading of the particle wave packet, 
and was derived by evaluating the disconnected part of the correlation function. 

Finally, regarding discussion of Sec. IIII Dl it should be stressed that for the purpose of experimental verification 
of Eq. (|3"3"1) only the genuine 1/f noise data was used, i.e., the data that fits the law /~ 7 in which 7 = 1, within 
experimental error. Otherwise the comparison would be meaningless, even for / « 1 Hz. Large deviations of 7 
from unity, observed in some thin films, are presumably due to back reaction of the conducting medium on the 
electromagnetic field produced by the charge carriers. This issue will be considered elsewhere. 
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APPENDIX A: GAUGE INDEPENDENCE OF THE CORRELATION FUNCTION 

Consider the theory of interacting scalar and electromagnetic fields described by the action 

S = Scf, + Sa , 

where is given by Eq. ([5]), and 

s A = -\f d^xF^F^ + S gf , S gf = ^f d A x (d^) 2 , Ffj, u = d^Av - d v A^ . 
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For arbitrary constant parameter a, the gauge fixing term describes the generalized Lorentz gauge. Let us introduce 
the generating functional of Green functions 



Z[J, 77, if] = J dAd<f>d<f>* exp ji (s + J d 4 :r[J / M, i + r?*0 + r)<f)*] J j , (Al) 

where J, 77, rj* denote sources for the fields A,<f>*,<f>, respectively. Vanishing of Z under the gauge variation of the 
functional integral variables 

8A I1 = d^(x) , 5(f) = ie£(x)<f) , 5<j>* = — ie^(x)(j)*, 
with £(x) a small gauge function, leads to the Ward identity 

- W^(y)Z + + ie^iv)^ - ieviv)^ = . (A2) 

Since we are interested in the connected contribution to the correlation function, we rewrite this identity for the 
generating functional of connected Green functions, W — —ihiZ, 



-d SW + ierfUj) — - ierjjy) — 



Vfo) + -d M7T77T + ietf^)^-^ ~ ieriiv)-^ = . (A3) 



The consequence of this equation we need is obtained by functional differentiation with respect to 77, rj* , and twice 
with respect to J, with all the sources set equal to zero afterwards, 

□» 6*W . (4) 6*W 



'SJMS^yM^WzWW '5J a (x)5J,3(x')6r)(z)6r,*{y) 
-ieS^Hv - z) — = 

Fourier transform of this identity with respect to y reads 



a J t 



SJniyWaWSJpWSriWSTriz') 



•i > -t S 4 W 
1 ; ( 5J Q (x)(5J / 3(a;')(57](z)V(^) ' 

The argument of the Fourier transform is purposely denoted here by k\ to stress that the left hand side of this equation 
corresponds to the variation of the Green function we dealt with in Sec. IIIII under gauge variation of the external 
field. Indeed, the longitudinal part of the photon propagator in the generalized Lorentz gauge has the form 

DU k ) = ~^- (A5) 

Therefore, contraction with the factor k\kx^/a is equivalent to amputation of the photon propagator attached to the y 
vertex, followed by contraction of this vertex with fci M . Exactly the same result is obtained under a gauge variation of 
the external field coming into this vertex. The only difference with the Green function we considered in Sec. IIIII is that 
the external scalar lines in Eq. (|A4[) are the particle propagators. To promote them into particle amplitudes, according 
to the standard rules, Eq. (|A4|) is to be Fourier transformed with respect to the variables z, z', and then multiplied by 
a(q)a* (q')(m 2 — q 2 )(m 2 — q' 2 ), where the arguments q, q' of the Fourier transformations with respect to z, z' are to be 
taken eventually on the mass shell. But these operations give zero identically when applied to the right hand side of 
Eq. l|A4p . because each of the factors e~ lklZ , e~ lklZ makes the corresponding particle propagator nonsingular on the 
mass shell. For instance, the first term in Eq. (IA4[) gives rise to the contribution of the form (m 2 — q' 2 )D^(q' + hi) 
times terms nonsingular on the mass shell. For k\ ^ 0, the function D^(q' + k\) is also nonsingular at q' 2 = m 2 , and 
hence this contribution vanishes on the mass shell. 

Thus, the correlation function is invariant under the gauge transformations of the external field, which are part of 
the gauge freedom in the theory. The other part is related to the explicit dependence of the photon propagator on the 
choice of the gauge conditions used to fix the gauge invariance of the action. As is well known, it is the longitudinal 
part of the propagator that depends on the gauge, and the most general Lorentz-invariant form of this part is given 
by Eq. (| A5|) in which a is to be regarded as an arbitrary function of k 2 . It is not difficult to see that variations of 
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a(k 2 ) do not affect the observable quantities. Recall, first of all, that we are interested ultimately in the fluctuations 
of gauge-invariant quantities such as the electric field strength. The a-independence of these quantities is a direct 
consequence of their gauge invariance, because variations of a(k 2 ) give rise to terms that are pure gradients with 
respect to the spacetime arguments x,x', as is easily verified by substituting the expression (|A5|) in place of one or 
two photon propagators in Eq. f|13[> . Then, if the vector potential contribution to the field strength is negligible, as 
is the case in our nonrelativistic calculation (recall the condition \q\ <C m used throughout), the voltage correlation 
function can be found by integrating the correlation function for the field strength with respect to x,x' using the 
relation E = — V A . 

Thus, gauge- independence of our results expressed by Eqs. (|32f . ([33)1 is proved. 
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Figure captions 

Fig.l: Feynman diagrams representing connected part of correlation function. Wavy lines denote photon propaga- 
tors, solid lines massive particle, q and p are the particle 4-momentum and 4-momentum transfer, respectively. 

Fig. 2: Symbolic diagrammatic picture of the effect of particle collisions and external electric field (dashed line) on 
the particle wave function. 

Fig. 3: Feynman diagrams describing the first order external field correction to the particle propagator. 




FIG. 3: Kazakov, Quantum fluctuations of Coulomb potential as a source of flicker noise. The influence of heat bath 



